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Abstract 

We propose an exactly solvable model to describe the properties of atomic 
thickness hybrid ferromagnet-superconductor-ferromagnet (F/S/F) struc- 
tures. We show that the superconducting critical temperature is always higher 
for antiparallel orientation of the ferromagnetic moments. However at low 
temperature the superconducting gap occurs to be larger for parallel orien- 
tation of the ferromagnetic moments. This leads to a peculiar temperature 
dependence of the proximity effect in (F/S/F) structures. 

PACS : 74.50+r, 74.80.Dm. 



I. INTRODUCTION 

In recent years the superconductor- ferromagnet (S/F) hybrid structures attract a steadily 
growing interest. The actual progress in the preparation of S/F/S junctions permits to ob- 
serve the transition from normal 0-junction to the so-called 7r-junction with the change of 
the thickness of the ferromagnetic layer [1,2]. Such a behavior has been predicted in [3,4] 
and it is related to the oscillations of a superconducting order parameter in the ferromag- 
netic layer. Another manifestation of the spin-dependent transport in S/F hybrid structures 
is the dependence of the critical temperature of metallic F/S/F sandwiches on the mutual 
orientation of ferromagnetic moments of the outer layers. On the basis of Usadel equa- 
tions, applicable in the dirty limit, it has been demonstrated in [5-8], that the antiparallel 
orientation of the ferromagnetic moments is more favorable for superconductivity, i.e. it 
corresponds to the higher superconducting transition temperature and to the larger ampli- 
tude of the superconducting order parameter. Recently this effect has been observed on 
experiment [9]. Note also that for the first time the coupling between ferromagnets through 
a superconducting layer has been treated theoretically in [10] for the case of ferromagnetic 
insulators and observed on experiment in [11,12]. 

At the same time in [13,14], a microscopic multiterminal model for S-F hybrid structure 
at T = has been proposed, where the ferromagnet was described by different density of 
states for opposite orientations of the magnetic moment. Then, it has been deduced that 
the superconducting gap is less influenced by the ferromagnet ism for the case of parallel 
orientation of the F electrodes. Hence at T = 0, the parallel orientation of the ferromagnetic 
moments appears more favorable for superconductivity. 
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In the present work we analyze the properties of the microscopic model system comprising 
three coupled atomic-scale layers: one superconducting layer in between of two ferromag- 
netic layers, which can have parallel or antiparallel orientation of the magnetic moments. 
This model may be relevant for the description of layered superconductors like RuSr 2 GdCu 8 , 
where superconducting and magnetic layers alternate [15] as well as artificial S/F structures 
obtained by molecular beam epitaxy. Such a model can be solved exactly and so we may 
verify that the results for the critical temperature for a diffusive regime of electrons motion, 
obtained in [5-8], are also qualitatively applicable for a clean atomic layered system. We 
demonstrate that in all cases the antiparallel alignment of ferromagnetic moments corre- 
sponds to the higher transition temperature to the superconducting state. However at low 
temperature the situation may be different and as it follows from our analysis, the super- 
conducting gap for parallel alignments of ferromagnetic moments increases faster with the 
decrease of the temperature comparing with the antiparallel alignment. At some temper- 
ature T* the gaps for parallel and antiparallel alignment coincides and at T < T* the gap 
is higher for parallel alignment. So at T = the parallel alignment is more favorable for 
superconductivity while near T c the situation is inverse. This result is in accordance with 
[13,14], where the superconducting gap have been calculated at T = in the framework 
of quite different model of bulk superconductor in the contact with small ferromagnetic 
electrodes. 

II. SPIN-ORIENTATION-DEPENDENCE OF SUPERCONDUCTING 
TRANSITION TEMPERATURE 

Adopting the model of [16], we consider a three layers system with one superconduct- 
ing layer between two ferromagnetic layers, see Fig. 1. It is supposed that the coupling 
between layers is realized via the transfer integral t, which is relatively small (t < T c ), so 
the superconductivity can coexist with ferromagnetism in the adjacent layers. Introducing 
the notations and 77+ for electrons creation operators in F layers and tp£ in S layer, the 
Hamiltonian of the system can be written as 

H — Ho + H^p + + 

Ho = J2 [&(p)vt(p)<p*(p) + m^(p)Mp) + C(pH(p)vAp)] , (i) 

p,a 

p 

Htf, = t^ bPt(p)¥v(p) + <pZ(p)il>*(p)] , 

H ^ = tJ2 [^(p)vAp) + v+(p)Mp)} , 

p,a 

where H describes the free electrons motion in F layers with the spectra £$(p) and Q(p) 
respectively, and with the spin-independent spectrum in S layer. The BCS pairing in 
the middle S layer is treated in in the mean field approximation [17] and the coupling 
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between neighboring layers via the transfer integral t is described by and H^ v . Note 
that the electrons spectra in (1) are calculated from the Fermi energy. The superconducting 
order parameter satisfies the usual self-consistency equation 



d 2 p 



2 ' 



(2) 



where, A is the Cooper pairing constant which is assumed to be non zero in S layers only. 
Introducing in the usual way [17] the Green functions 



G a f3 



< T T (^+) > , F+ =< T T (^+) > , 

< T T (^J) > , F$ =< T T (ifcty+) > , . 



and writing the corresponding equations for the Green functions, we finally obtain the 
following exact expressions for the Green functions in S layer: 



F + - 



ab* + | A 
A* 



ab* + I A 



2 ' 



2 ' 



G 
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F + 



a*b+\A 
A* 



2 ' 



(3) 



a*b + |A 



2 ' 



where 




b = iuo - £ - t 2 

The self-consistency equation for superconducting order parameter is now written as 

d 2 p 



A* = |A|T^ j 



F + - 



(2tt)' 



To calculate the critical temperature of the superconducting transition it is sufficient to know 
the anomalous Green function F^ in the linear approximation on A*. Then it is convenient 
to write the linearized self-consistency equation for T c in the following form [17] 



-©-•Sin 



Mf) # 



(4) 



where we have defined a reduced function / = 1/ (a*b) and T c o is the bare mean-field critical 
temperature of the central S layer in the absence of the proximity effect (i.e. for t — 0). 

In the case when both ferromagnetic layers are equivalent we may distinguish two differ- 
ent situations 

either the parallel (P) orientation of the magnetic moments, where £j = £j? = £j? and 

£i = £|. — Cj. > 
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or the anti-parallel (AP) orientation, where £|- = £j? = £j and £j_ = ^ — ^. 
Firstly we consider the situation when the electron dispersion spectra in the ferromagnet 
differs only by the exchange field h from the electron spectrum in S layer i.e. £j = £ — h and 

£j C ■ I' 

Performing in (4) the integration over energy £ and summation over Matsubara's fre- 
quencies uj = {2k + 1)"kT c , we finally obtain 



In 



T, 



cO 



2t 2 



T cP n 2 T 2 8t 2 + h 2 



h 2 ^ 



h 



2nT, 



cP 



+ 



(5) 



At 2 



- , V8t 2 + h 2 + h\ ^ ( V8t 2 + h 2 -h 

$1 7T-^ + $1 



2ixT cP 



2ixT cP 



In 



cO 



2t 2 



T cAP n 2 T 2 2t 2 + h 2 



h 2 ^ 



V2t 2 + h 2 
2nT cAP 



+ 2r$i 



VWTh 2 

kT cAP 



where T cP (T cAP ) is the superconducting transition temperature for parallel (antiparallel) 
orientation of the F-layers magnetic moments and we define the function <£>! through the 
Digamma function ^ as (7 is Euler constant) 



$i(x) = ^| 7 + 21n(2) + i 



IX 



IX 



Using (5) we can determine the relative variation of critical temperature 5T = (T c — T c0 )/T c0 
in the two followin g limi ts : 

i) in the limit — > up to the fourth order over t we have: 



5T P ~ 



7C(3) 



(2vrT c0 ) 
1 



+ 



(27rT c0 ) t 



(2ttT c0 
1085 



62C(5) - ^C(3) 2 ) t 4 + ^C(5)*V 



+ 



lg-C(3)C(5) - ^C(7)j t 4 /* 2 - ^C(7)t 2 /. 4 



1397C(7) t 4 /i 2 
512tt6 if" 



5T P - 5T AP ~ 



t 4 /i 2 



-0.0029- 



T 2 ' 



(6) 



The last result, the difference STp — 5T A p, was already found for S/F multilayer in [7] (after 
correcting a sign mistake). 

ii) in the limit ^ t2 + h2 — > 00 and for t « h we obtain : 



5T P ~ -4 



7? 



7 + In 



cO 



+ 



7C(3) t 2 
4^ 2 T 2 



5T P - 5T AP0 



7C(3) t 4 
4vr 2 h 2 T 2 



~ -0.21 



(7) 



Also the difference of the critical temperatures 5T P — 5Tap coincides with that of the S/F 
multilayer [7]. 

We see that in all cases the difference between critical temperatures for parallel and 
antiparallel alignment is proportional to t 4 , while the decrease of the critical temperature 
itself is proportional to t 2 . This fact demonstrates that the spin-orientation dependence 
of T c is related with a rather subtle interference effect between electrons coming from 
ferromagnetic layers. 

Now we consider the case of a ferromagnetic half-metal, which we may model in (3) 
by £j. = £ + h and £j — > oo, that corresponds to a zero density of sates for the electrons 
with spin orientation along magnetic moment in ferromagnetic layers. Note that the limit 
£l - ► oo is equivalent to take the transfer integral for electrons with spin up orientation 
equal to zero. This may be simply demonstrated from the initial equations for the Green 
functions. Performing integration over £ and summation over Matsubara frequencies in the 
self-consistency equation, we can deduce the relative temperature variation 5T = (T c — 
T c0 )/T c0 at the order up to t 4 , in terms of the dimensionless field h = h/(2nT c0 ) 



ST P = -Q 1 [h 



+ 



7T 2 T 2 
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vr 4 T 4 ' 
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where the function <&i was defined before and the function $ 2 is 



$J X ) = — 
2V ; 16i 
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In the limiting case h 
have 



0, which corresponds to the situation h -C T c0 ,or simply £j = £, we 



5T P = - 



:C(3) 



vr 2 T c 2 



+ I ||C(5) - I (|C(3) 



and the critical temperature difference is 
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8T P -8T AP = -_ C (5)-^ ^ -0.50 
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In the opposite limit /i ^> T c0 , we find 

1 



5Tp 
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and 



ST P - ST AP i^ln ( -A_ ) . (10) 
h 4 \2nT r 1 



-cO 



Then we may conclude that in all cases the critical temperature is higher for the antiparallel 
alignment of the ferromagnetic moments. 

III. SUPERCONDUCTING GAP AT LOW TEMPERATURES 



In this section we demonstrate that the proximity effect at low temperature is very special 
and that the superconducting gap at T = 0, for half-metal and for usual ferromagnet with 
small exchange field, is higher for parallel orientation in accordance with [13,14]. 

Firstly we consider the model of a half-metal with £j — > oo and ^ = £. So the self 
consistency equations for parallel and antiparallel orientations may be written as 



d£ 



(11) 



where iV(0) is the electron density of state. 

In the limit of weak interlayer coupling t « T c0 and at T = 0, we obtain the following 
expression for the superconducting gap for parallel orientation Ap 



2tt 



t 2 

1 + 2- y exp(2^) 



dO, (12) 



where Ao is the gap of an isolated S layer. Performing expansion over t, it may be demon- 
strated that at all order over t 2 the corrections to Ap disappear, and Ap = A . So the 
proximity effect for the case of half-metal vanishes. It may be understood as an impossibil- 
ity of Cooper pair destruction. Indeed at T = 0, the disappearance of the Cooper pair in 
S layer means that two electrons with opposite spins must leave it. Due to the insulating 
character of neighboring F layer for some one spin orientation it becomes impossible and so, 
the Cooper pair is not destroyed at all. 

On the other hand, for the antiparallel alignment this argument does not works and the 
proximity effect leads to a decrease of the gap Aap 

Therefore the superconducting gap variation is (A — Aap) /A ~ 2 (t/A ) ln(A /t) . As 
the superconducting transition temperature is higher for the AP case, it means that above 
some temperature T* the gap for parallel orientation becomes smaller than that for antipar- 
allel orientation. 
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To find this temperature T* in the limit t « T c q we perform an expansion over t in the 
self-consistency equations for A^p = Ap = A*. After integration over energy £ we have the 
following equation for the ratio X = A*/ (nT*) 

2{2k + l) 3 + 5(2k + 1)X 2 - 2[5(2A; + l) 2 - X 2 }[(2k + l) 2 + X 2 ] 1/2 
V (2fc + l)[(2fc + l) 2 + X 2 ]3/2[(2A; + 1) + ^(2^ + l) 2 + X 2 ] 4 " ' (M) 

from what we find A*/T* = 4.22. Taking into account that in the limit t << T c0 the 
temperature dependence of the superconducting gap in the first approximation is the same 
as that for an isolated S layer A (T), we directly find the temperature of gap inversion 
T* = 0.41T c0 . 

Similar calculations can be performed for the standard ferromagnetic case with £j = £ — h 
and £i = £ + h in the limit t,h <C T c0 . In the result we obtain the following equation for the 
ratio X = A*/(nT*) 



E 



f -2{2k + l) 4 - 16(2fc + 1) 2 X 2 - 6X 4 

\ (2fc + 1) 2 [(2A; + l) 2 + X 2 ] 3 / 2 [(2A; + 1) + ^{2k + l) 2 +W] 

+38{2k + l) 4 - U(2k + 1) 2 X 2 - X 4 



= 0, (15) 



(2k + l) 3 [(2k + l) 2 + X 2 ][(2A; + 1) + ^{2k + l) 2 + X 2 ] f 
which gives A*/T* = 3.61, and the temperature T* of the gap inversion is T* = 0.47T c0 . 



IV. CONCLUSION 

In the present work we have demonstrated that, in the framework of the considered 
microscopic model, in all cases the superconducting transition temperature is higher for 
antiparallel orientation of the ferromagnetic moments. This is consistent with the predictions 
made in [5-8], on the basis of the diffusive dirty limit model and recent experimental results 
[9]. We may consider this result as a quite general one and model independent. On the other 
hand the superconducting gap at low temperature for the case of half-metal and for usual 
ferromagnet with small exchange field occurs to be larger for parallel orientation. This is 
in opposite with the diffusive model prediction [8], but in accordance with the T = result 
[13,14] obtained in the framework of the multiterminal model for S-F hybrid structures. It 
may be interesting to calculate the superconducting critical temperature in the model [13,14] 
to compare with our predictions. 

Finally we conclude that the spin-orientation dependence of the superconductivity in 
F/S/F atomic layers structure may be quite different at high and low temperature regimes. 
Note also that the condition of small interlayer coupling t « T c q, may be crucial for the 
conclusion that there is a gap inversion for parallel and antiparallel orientation with tem- 
perature decrease. Indeed, it has been demonstrated in [18], that the strength of interlayer 
coupling may qualitatively change the critical temperature dependence via the exchange 
field. 

Note also, that strictly speaking, the superconductivity is impossible in atomic 2D sys- 
tem. However, our results must be qualitatively applicable for systems consisting of several 
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consecutive S and F layers, as well as for S/F multilayered systems, where the fluctuations 
are strongly suppressed. 
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Figure caption 

FIG. 1. Geometry of three layers F-S-F system, one superconducting layer is between 
two ferromagnetic layers. 
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Figure 1 



